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Abstract: We deal with scalar field coupled to gravity in five dimensions in warped 
geometry. We investigate models described by potentials that drive the system to support 
thick brane solutions that engender internal structure. We find analytical expressions for 
the brane solutions, and we show that they are all linearly stable. 
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Branes in higher dimensional theories provide an interesting new procedure for resolv- 
ing problems involving cosmological constant and hierarchy [Q, |2| . The brane scenario that 
we are interested in the present letter is driven by a single real scalar field. It engenders 
the interesting feature of localizing massless graviton in the brane, efficiently reproducing 
four dimensional gravity in the brane. The study of scalar fields coupled to gravity in 
warped geometries has gained renewed attention recently ^, |6|, ^. In particular, in 
Ref. one investigates the splitting of thick branes due to a first-order phase transition 
in a warped geometry. The work of Ref. considers a complex scalar field coupled to 
gravity, and it shows that when the temperature approaches its critical value, an interface 
that interpolates two bulk phases breaks into two separated interfaces, giving rise to a new 
phase which appears in between these two interfaces, leading to an effect which is known in 



condensed matter as complete wetting — see also Ref. |10] for other details in the subject. 

In the present work we show that similar results also appear at zero temperature, in 
a simpler model, which depends on a single real scalar field, coupled to gravity in five 
dimensions in a warped geometry. The results seem to be of compelling interest to high 
energy physics, and are depicted in the figures below. The potential that controls the scalar 



field depends on odd integers in a very specific manner |11], allowing the appearance of 
thick branes that host internal structure, in the form of a layer of a new phase squeezed in 
between two separate interfaces. The potential that we investigate here was first studied 



in Ref. |11|, and there it was shown that although it depends on a single real scalar field, it 
supports 2-kink solutions of the BPS type, which are stable solutions that solve first order 
differential equations and engender the feature of being composed of two standard kinks. 
These topological defects are richer than the standard topological kinks, and they may be 
seen as defects that host internal structures. This feature establishes the main motivation 
of this work, in which we explore the possibility to offer a brane scenario where the brane 
entraps internal structure. The idea is similar to the investigations done in |p!^, 13 1 in ffat 



spacetime, in models described by two real scalar fields, and in in curved spacetimes, 
in brane scenarios involving two higher dimensions. In the present work, however, we shall 
deal with very specific models, described by a single real scalar field, which couples with 
gravity in warped spacetime in one higher dimension. 

To show how the above potential behaves in the brane scenario, we follow Refs. @, g 
and we consider the action 



S = J d X dy^/\( 

where g = det{gab) and the metric is 

ds"^ = Qahdx^dx^ = e2^r/^.^dx^(ix^ - dy'^ (2) 

where a, 6 = 0, 1, 2, 3, 4, /i, = 0, 1, 2, 3, 'q^y = (1, —1, —1, —1) and e^"^ is the warp factor. 
We suppose that the scalar field and the warp factor only depend on the extra coordinate 
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y. In this case the equations of motion are 

," + 4^V^| (3) 
A" = -1<P" (4) 
A" = - lv{^) (5) 
where prime stands for derivative with respect to y. Our model is given by 



where [11 1 



, , 2v 2p-i 2n 2p+i 

WM) = ^,^-^-^,*-^ (7) 

The parameter p is odd integer and controls the way the scalar field self-interacts. For 
clarity, in Fig. [1] we depict the potential Vp{(j)) for p = 1,3, 5. 

The function VFp(</>) was first introduced in Ref. |11| in flat spacetime. It was inspired 



in the work on deformed defects |15|, and it appears with the help of the deformation 



prescription there proposed. To see how this is done, we suppose that the spacetime is flat. 
In this case the potential is given by V{(j)) = {l/2){dW /dcj))'^ . We consider the function 



W{(j)) = (j) — (j) /3, which identifies the standard (j) model. We now follow |15| and introduce 



the function f{4>) = (^^^^ , which is well-defined for every 0, for p odd integer. This function 
can be used to change VF(</>) = (t)—(p^/3 to a new model, defined by Wp{(j)). The deformation 
prescription leads to the general result (df /d(j)){dWp/d(j)) = {dW /d(j))[(j) — > where in 

the right hand side one needs to change (p by after taking the derivative of VF((/>). 

For the case under consideration we get Wp{(j)) = {p/2)Wp{(l)), that is, {p/2) times the 
expression in Eq. (0), which is used to define the potential Vp{4>). 

We return to curved spacetime, and we recall that potentials defined as in Eq. (|6|), 
through the introduction of the function W = W {</)), appear very naturally in supergravity, 
and there W is named superpotential — see Refs.[16|, 17, 18| for more details on this subject. 

The use of Wp{(j)) to define the potential provides an important approach, which leads 



to first-order equations that solve the equations of motion — see for instance Refs. [19, 20 
|6l. The first-order equations are given by 

^' = i^^ (8) 

^ 2 d(P ^ ' 

and 

A' = --Wp (9) 

The first-order equation for cp reproduces the first-order equation for the scalar field in flat 
spacetime [^], apart from the numerical factor of 1/2. We solve this equation to get |11| 



.p{y) = tanhP { ^ ) (10) 
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Figure 1: Plots of the potential Vp{(j)) for p — 1,3, and 5. Here and in the other figures the 
thickness of the lines increases with increasing p. 



We see that for p = 1 we get the standard solution. However, for p = 3,5,... we get 



2-kink solutions, which we show in Fig. [2] — see Ref. [11| for more details on such defect 
structures. 




Figure 2: Plots of the solutions (j>p{y) for p — 1,3, and 5. 



The first-order equation (|9|) for the warp factor is new for Wp given by Eq. (|^, and 



we have been able to solve it for Ap{y) analytically. We get 



My)-- 



1 p 



3 2p+ 1 
In 



tanh^P 



cosh ( - 



_- 
p) 3 

y-tanh2" (y 

n=l ^ 



2p-l 2p+l 



We plot the warp factor e^'^p^^^ corresponding to these solutions for p = 1, 3, 5 in Fig. [3]. 

The scalar matter field and the warp factor contribute to the brane scenario according 
to the solutions (|^) and (|ll|). We use these results to obtain the matter energy density 
^00 (y) ^ ^ function of the extra dimension. The analytic expression for Tqq is rather 
awkward, thus we decided to plot it in Fig. [4] to show the matter behavior when 0(y) and 
A{y) solve the first order equations (S) and (0), respectively. 




Figure 3: Plots of the warp factor e^^'''^^^ for p = 1,3, and 5. 

The plots depicted in Fig. [4] show that the matter field gives rise to thick branes 
composed of a single {p = 1) or two {p = 3,5, ...) interfaces. In the last case, for p = 3, 5, 
one sees the appearance of a new phase in between the two interfaces where the energy 
density of the matter field gets more concentrated. The scenario for p = 3,5, is that of 

We use Figs. [3] and [4] to see that 



a brane that supports internal structure [12| , 
the warp factor is essencially constant (for p 



= 3,5,. 



in the very inside of the brane, in 
the region where the energy density vanishes. This unveils the presence of a new phase 
inside the brane, protected by the two interfaces that appear centered at the points which 
maximize the energy density. For p = 1 the scenario is quite different, and the brane does 
not support the internal structure that we clearly see for p = 3,5, ... 

The scalar matter model (0) that we have used leads to the solutions ([lO|) and (pi]), 
which are depicted in Figs. [2] and [3], respectively. We go further into the subject and we 
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Figure 4: Plots of the matter energy density TgQ(i/)/10, Toq(i/), and T^Q^y) for p = 1,3, and 5. We 
have plotted 1/10 of Tqq to ease comparison. 



examine stability of the above braneworld scenario. The key issue here is stability of the 
geometry induced by the above solutions. We do this perturbing the metric, using 

ds^ = e^^^y^ [t]^, + h^,)dx^'dx'' - dy^ (12) 



and supposing that h^^ = h^i^{x,y) represent small perturbations. Here we follow |^] and 
we choose hf^i, in the form of transverse and traceless contributions, h^u. In this case one 
gets the equation 

h';,, + AA'h'^^ = e-^^ah,, (13) 

which describes linearized gravity. Here □ stands for the wave operator in (3, 1) dimensions. 
We notice that in this sector, gravity decouples from the matter field. 

We use the new coordinate z, which is defined by dz = e~^^y^dy. Also, we set 

V(x,z) =e*^-^e"i^(^)i?^,(z) (14) 

In this case the above equation (|l3l) becomes the Schrodinger-like equation 

dz^ 

where the potential is given by 



+Up{z)H^, = k^H^, (15) 



U,{z) = ^A;{z) + ^-A'^\z) (16) 



The Hamiltonian in eq. (^) can be written in the form 

d_ 3 ,\ f d_ 3 

and this ensures that k is real, and /c^ > 0; thus, there is no unstable tachyonic excitation 
in the system. We solve for the zero modes {k = 0) to get 

HlP){z)=N(^Je'^^(^y' (18) 
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where Nj^J is a normalization factor. In Fig. [5] wc plot the potentials [/i(z)/10, ^73(2;), and 
U^{z), and the corresponding zero modes for p = 1, 3, 5. In Fig. [6] we plot wave functions 
for massive modes in the two regions where fc^ < max[C/p(2;)] and k"^ > m.SiX.[Up(z)]. These 
wave functions confirm the expected result, that they describe motion in the bulk, not 
bound to the brane. 




Figure 5: Plots of the potentials Ui {z)/ 10, U3{z), andU^lz) (left) and the corresponding zero 
modes (right). We have plotted f7i(2;)/10 to ease comparison. 




Figure 6: Plots of wave functions for nonzero modes for < ni&x.[Up{z)] (left) and for k^ > 

max[[/p(z)] (right) for p = 1,'S, and 5. 



In the brane scenario that we have just examined, in order for the zero modes to 
describe localized four dimensional gravity, normalizability is essential. To ensure normal- 
izability, the zero modes as a function of z must fall off faster than For this reason, 

in Fig. [7a] we have plotted the zero modes and their intersection points with the function 
f{z) = z~^/^. And in Fig. [7b] we show that the intersection points Zp increases with p 
with decreasing derivative. This result indicates that the zero modes always fall off faster 
than showing that they are normalizable, giving rise to localized four dimensional 

gravity for every p odd integer. 

A consequence of this result is that the massive modes are strongly suppressed. The 
structure of the volcano potentials plotted in Fig. [5a] could contribute to the appearance 
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of resonances for k < max[C/p(z)]. To check for the presence of resonances, we have inves- 
tigated massive states for several different values of k in the interval [0, max[?7p(2;)]], for 
p = 3 and 5. We have found no evidence for resonances. Thus, our model seems to support 
no resonance for p odd integer, despite the fact that the values p = 3,5,... lead to the 
generation of two interfaces which inhabit the thick brane presented in this work. 
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Figure 7: Plots of the zero modes for p — 1,3,5 (left), and intersection points with f{z) = z^^/^, 
which is shown as the dashed line. We also plot (right) the intersection points between the zero 
modes and the function /(z) = z^^^^ for several values of p. 



Let us now comment a little further on the 2-kink solutions that we have found in 
Eq. ([lO|). We can have a better view of these defect structures calculating y^, the position 
which identifies the center of the defect. We identify yp with the maxima of the derivative of 
the defect solutions, that is, yp = max[d(f)p{y) / dy]. We get yp = itparctanh \/ {p — l)/{p + 1). 
For p = 1 it gives yi = 0, as expected, but for p = 3,5, we obtain two non-zero values, 
which are directly related to the center of the interfaces that inhabits the defect. In flat 
spacetime, the values yp exactly identify the maxima of the energy density 

1 /d(hr,V 1 fldW, 



for the corresponding solutions (pp{y) — the unusual 1/2 factor that multiplies Wp appears 
because in flat spacetime we have to redefine Wp as Wp/2, in order not to change the first- 
order equation (|8|) for the scalar field (j). This result shows that the 2-kink solutions are 
composed of two interfaces, each one being centered at each one of the two points yp that 
we have just obtained. In curved spacetime, however, the position of the maxima of the 
matter energy density is different. To see this, we use TQQ{y) to identify the maxima of the 
matter energy density in curved spacetime. The function T^Q{y) is depicted in Fig. [4] for 
some values of p. There are maxima at the points y^ = 2.42, y^ = 5.36, yj = 8.70, which 
are different from the points y^ = 2.64, y^ = 5.73, y^ = 9.22, which identify the maxima 
of the energy density in fiat spacetime, or the center of the interfaces of the 2-kink solutions 
that we have already found. This result leads to the conclusion that the curved spacetime 
induces an attraction between the interfaces that appear from the 2-kink solutions. We 
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have investigated this effect for several p, and we could verify that the relative attraction 
decreases for increasing p. 

We recall that solutions similar to the 2-kink defects depicted in Fig. [2] have already 
been considered before in flat spacetime. See for instance Ref. which investigates a 
model described by a complex scalar field. We can also deal with two real scalar fields 



and consider the model investigated in |21| in flat spacetime. It is defined by VF(</), x) = 
(j) — 4>^ /?> — rcjyx^, where r is a real parameter. It supports kink-like solutions p2| , [2^ of the 
form presented in Fig. [2]. This model was very recently used to describe matter coupled 



to gravity in warped geometry |24], with the result that the warp factor resembles the form 
depicted in Fig. [3]. The models investigated in the present work are similar but simpler, 
since they are described by a single real scalar field. 



The model investigated in |21| support two-field solutions in the sector defined by the 



minima (±1,0) in the {(j), x) plane |25|. We may couple this model to gravity, to investigate 
whether these two-field solutions generate braneworld scenarios. Other extensions may 
involve scalar fields described under the general guidance of the Z2 x symmetry — see 
for instance the second work in Ref. |25| — and they may also lead to issues involving 
branes that host network of non trivial structure. 

We would like to thank F.A. Brito and J.R. Nascimento for helpful discussions. We 
also thank PROG AD and PRONEX for financial support. DB and CF thank CNPq for 
partial support and ARG thanks FAPEMA for a fellowship. 
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